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FIG. 1. Average bulk modulus B{e,L) of randomly site- 
displaced Penrose approximants with period L and disorder 
strength e, under Periodic (PBC, a) and b)) and Fixed (FBC, 
c) and d)) boundary conditions. Lines in a) and c) are guides 
to the eye. The thin line in b) is i? ~ Thick lines in b) 

and d) are, respectively, B ~ {e^ and B ~ 


The discovery of isostaticity in sphere packings [I| and 
network glasses UM has inspired a great deal of activity 
in the field of isostatic networks. Recent studies Si sug¬ 
gest that all elastic modulii of geometrically disordered 
isostatic networks go to zero with increasing linear size L, 
if disorder is uncorrelated. Packings of hard frictionless 
spheres or discs, on the other hand, have nonzero com¬ 
pressive modulus B ii , despite being isostatic i and 
disordered, because their contact network is not random, 
but tuned to avoid negative forces. Contact disorder is 
correlated in these systems. Attempts to model sphere 
packings as randomly disordered isostatic networks have 
therefore failed. However, in a recent letter Q, Stenull 
and Lubensky (SL) claim that randomly disordered Pen¬ 
rose networks have nonzero B for large L. The present 
numerical study, using high-precision Conjugate Gradi¬ 
ent to solve the elastic equations shows their claim to 
be incorrect, and clarifies the reason for their misinter¬ 
pretation of the data. Figs. [1^ and b show B{e,L) for 
Penrose periodic approximants of orders 5 to 12 (up to 
8 X 10^ sites), whose sites are randomly displaced within 
a circle of radius e. B behaves roughly as 1/L^ for large 


L (Fig. [TJi). However, because B{e = 0,L) = OVL B 
grows as Q when << 10^. The asymptotic 
regime L » Lq m (10/e)^/^ is hard to reach for small 
e. This has been noted already @ for other disordered 
isostatic netwoks. The data reported by SL @ (derived 
from normal-mode calculations for a single, unspecified, 
value of e) are similar to our results for e = 10“^ in 
Fig. [T^, i.e. B appears to saturate. Our scaling analy¬ 
sis in Fig. [TJ) shows that this is a finite-size effect: the 
true asymptotic behavior B L~^ would only be seen 
at much larger sizes for this value of e. Further valida¬ 
tion of our claim that B ^ 0 for large L is provided by 
the following: fixing a line and a row of sites produces 
Penrose networks with Fixed BC. B^^'^ is seen to go to 
zero with size when e^L » 1 (see Figs. [TJ: and d). But 
B^bf jg ^ rigorous upper bound for B^^^. Therefore, 
j^PBC Q Jqj. ig^]-gg i, ag -vvell. We additionally mention 
that the effects of geometric disorder on elastic constants 
can be described analytically for small e, giving rise to ra¬ 
tional expressions for B(e, L), that predict an asymptotic 
power-law behavior B{e,L) ^ L~^ when e ^ 0. Details 
will be provided somewhere else 0 . We conclude that the 
main point raised by SL is not justified: generic Pen¬ 
rose networks with uncorrelated geometric disorder have 
zero bulk modulus for large sizes. They are, therefore, 
no better suited to model jammed packings than any of 
the previously studied isostatic networks with geometric 
disorder. The authors thank CGSTIC of CINVESTAV 
for computer time on cluster Xiuhcoatl. 
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